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Using recently developed quantum SU(2)x U(l) rotor approach, that provides a self-consistent 
treatment of the antiferromagnetic state we have performed electronic spectral function calculations 
for the Hubbard model on the square lattice. The collective variables for charge and spin are 
isolated in the form of the space-time fluctuating U(l) phase field and rotating spin quantization 
axis governed by the SU(2) symmetry, respectively. As a result interacting electrons appear as 
composite objects consisting of bare fermions with attached U(l) and SU(2) gauge fields. This 
allows us to write the fermion Green's function in the space-time domain as the product CP 
propagator resulting from the SU(2) gauge fields, U(l) phase propagator and the pseudo-fermion 
correlation function. As a result the problem of calculating the spectral line shapes now becomes one 
of performing the convolution of spin, charge and pseudo-fermion Green's functions. The collective 
spin and charge fluctuations are governed by the effective actions that are derived from the Hubbard 
model for any value of the Coulomb interaction. The emergence of a sharp peak in the electron 
spectral function in the antiferromagnetic state indicates the decay of the electron into separate spin 
and charge carrying particle excitations. 



I. INTRODUCTION 

Recent high-resolution angle-resolved photoemission 
spectroscopy (ARPES) studies revealed a complicated 
character of electronic structure and quasiparticle (QP) 
spectra in copper oxide superconductors^ In the discus- 
sion of photoemission on solids, and in particular on the 
correlated electron systems, the most powerful and com- 
monly used approach is based on the Green's— function 
formalism. In this context, the propagation of a sin- 
gle electron in a many-body system is described by the 
time-ordered one-electron Green's propagator. A com- 
mon approach in describing strong electron correlations 
is based on consideration of the Hubbard modeL 2 It al- 
lows to study a moderate correlation limit observed ex- 
perimentally in cuprates and more consistently takes into 
account subtleties of the electronic structure, in particu- 
lar, a spectral weight transfer^ In qualitative sense the 
Hubbard model serves as the standard model of corre- 
lated electron systems, and has the same conceptional 
importance for interacting electrons as the Ising model 
for classical statistical mechanics. Here, the low num- 
ber of explicit parameters provides the ideal condition 
for a thorough test of the power and quality of analyt- 
ical and numerical methods. As a matter of fact in- 
tensive studies on this model have revealed subtlety of 
the results and controversies depending on approaches, 
and approximations. The two-dimensional (2D) Hub- 
bard model was studied in Ref. 0. The single particle 
self-energy was calculated by perturbation theory with 
the strength of the local Coulomb interaction as the ex- 
pansion parameter to the second order. The main effect 
of correlations is the transfer of the spectral weight to 
the high energies. The 2D Hubbard model on the square 
lattice was also studied in the presence of lattice dis- 
tortions in the adiabatic approximation. 5 In the absence 



of distortions the weight of the logarithmic singularity 
which characterizes the free system is reduced. Large U 
values give rise to the typical two-peaks situation corre- 
sponding to the infinite U limit. In other works, the 2D 
Hubbard model was considered with the quantum Monte 
Carlo method^^ recently also in the dynamical clus- 
ter approach^ To solve the cluster problem the Hirsch- 
Fye quantum Monte Carlo method has been combined 
with the maximum entropy method to calculate the real 
frequency spectra. The effect of larger clusters and in- 
teractions was also explored however, these results were 
restricted by computational limitations including, espe- 
cially, the minus sign problem. 10 Method based on nu- 
merical simulations for finite clusters precludes, however, 
to study subtle features of QP spectra due to poor energy 
and wave- vector resolutions in small size clusters* 1 ^ Thus 
careful analyzes of finite-size effects are important in nu- 
merical studies, especially for low-energy excitations. In 
the dynamical mean field theory (DMFT) the self-energy 
is treated in the single-site approximation which is un- 
able to describe wave-vector dependent phenomena^ 2 - 
To overcome the difficulties present in DMFT, various 
types of the dynamical cluster theory were developed. 13 
In these methods only a restricted wave-vector and en- 
ergy resolutions can be achieved, depending on the size 
of the clusters, while the physical interpretation of the 
origin of an anomalous electronic structure in numeri- 
cal methods is not straightforward. In the Two-Particle 
Self-Consistent (TPSC) approach ; 14 ' 15 that is based on 
enforcing sum rules and conservation laws, rather than on 
diagrammatic perturbative methods, a Luttinger-Ward 
functional is parametrized by two irreducible vertices 
that are local in space-time. This generates random 
phase approximation-like (RPA) equations for spin and 
charge fluctuations. The approach has the simple phys- 
ical appeal of RPA but it satisfies key constraints that 
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are always violated by RPA, namely the Mermin- Wagner 
theorem and the Pauli principle. The inherent difficulty 
of dealing with Hamiltonians appropriate for strongly 
correlated electronic systems originates from the non- 
perturbative nature of the problem and the presence of 
several competing physical mechanisms. In the fermionic 
systems with spin and charge excitations, however, the 
situation is much more complicated because dynamic 
quantum fluctuations are important even in large dimen- 
sions. Moreover, short-range spatial correlations play a 
key role, in particular magnetic correlations, leading to a 
strong tendency towards the formation of singlet bonds, 
as well as pair correlations. These correlations deeply af- 
fect the nature of quasiparticles. A new approach may 
therefore be needed, and a necessary requirement for a 
theory aiming to capture the essential physics of elec- 
trons correlated systems might be the inclusion of the 
fundamental ingredients of the physics as the underly- 
ing spin and charge symmetries as well as the associated 
ordered states. In this context the spectral properties of 
the two dimensional systems with strong magnetic fluctu- 
ations were investigated within the spin-fermion models 
in the quasistatic approach that neglects the effect of dy- 
namic spin fluctuations. The latter approach yields the 
two-peak structure of the spectral function for the anti- 
ferromagnetic correlations. 

The other completely different approaches rest on 
gauge theories that arise in models of strongly interact- 
ing electrons: it is often convenient to change variables 
from electron operators to other degrees of freedom that 
represents the intrinsic symmetries of the system under 
study. This route turns out to be beneficial since sev- 
eral quantum phases of matter (for instance a antifer- 
romagnet or a superconductor) may be characterized in 
terms of the symmetries that are broken spontaneously 
in that state . 17 ' 18 We employ this idea in the present 
paper to perform the one-particle spectral function cal- 
culations for the Hubbard model on the square lattice 
using a recently developed the quantum SU(2)xU(l) ro- 
tor approach^ The collective variables for charge and 
spin are isolated in the form of the space-time fluctuat- 
ing U(l) phase field and rotating spin quantization axis 
governed by the SU(2) symmetry, respectively. As a re- 
sult interacting electrons appear as a composite objects 
consisting of bare fermions with attached U(l) and SU(2) 
gauge fields. This allows us to write the fermion Green's 
function in the space-time domain as the product of the 
complex-projective (CP 1 ) propagator (which results from 
the SU(2) gauge fields), U(l) phase propagator and the 
pseudo-fermion correlation function. The problem of cal- 
culating the spectral line shapes now becomes one of 
calculating the convolution of spin, charge and pseudo- 
fermion Green's functions. The collective spin and charge 
fluctuations are governed by the effective actions that are 
derived from the Hubbard model for any value of the 
Coulomb interaction. We show that the method is use- 
ful for explicit calculation of spectral properties, enabling 
systematic inclusion of fluctuations of the charge as well 



as of the variable spin quantization axis. We demonstrate 
that the emergence of a sharp peak in the electron spec- 
tral function in the antiferromagnetic state points to the 
electron decaying into separate spin and charge carrying 
particle excitations. 

The paper is organized as follows. After introduction 
of the model Hamiltonian in Section II, we present in 
Sec. Ill the transformations to the phase and spin an- 
gular variables that reflect the basic symmetries of the 
Hubbard model. Section IV is devoted to the deriva- 
tion of the effective actions that govern the behavior of 
the system in the charge, spin and pseudo-fermion sec- 
tors, respectively. The self-consistent equations for the 
effective parameters of the model that follows from the 
effective actions are summarized in Section V. Section VI 
is devoted to a detailed analysis of the electronic spec- 
tral functions. Conclusions and discussions are given in 
Section VII. A number of technical details that pertain 
to the derivation of spectral functions is relegated to the 
Appendices. 

II. THE MODEL 

A prototype of theoretical understanding of the physics 
of strongly correlated systems is achieved by using simpli- 
fied lattice fermionic systems, in particular, the Hubbard 
model given by the Hamiltonian H = TCt + TCu- 

K = -tY, [4Wc Q (r')+h.c.]- M ^n(r)+^, (1) 

(rr'),a r 

where the Hubbard interaction term is given by 

Hu = U ni(r)H(r) (2) 

r 

and n(r) = n.f(r) + nj.(r) is the number operator. Here, 
(r,r') runs over the nearest-neighbor (n.n.) sites, t is 
the hopping amplitude, U stands for the Coulomb repul- 
sion, while the operator c' a (r) creates an electron with 
spin a =f (= 1),J, (= 2) at the lattice site r, where 
n a (r) = cJ l (r)c Q (r). The chemical potential fi controls 
the average number of electrons. The kinetic energy op- 
erator Ht reflects the electrons' itinerant features while 
the interaction operator Tiu forces the electrons' corre- 
lated motion or even their localization. 

Since the partition function often serves as a starting 
point for the calculation of thermodynamic properties, it 
is instructive to take a closer look at how this quantity 
may be obtained within the path integral formalism. To 
this end it is customary to introduce Grassmann fields, 
c a (rr) depending on the "imaginary time" < r < fj = 
1/ksT, (with T being the temperature) that satisfy the 
anti-periodic condition c a (vr) = — c Q (rr + (3), to write 
the path integral for the statistical sum 

Z = f [VcDc] e~ s[ ^ c] (3) 



3 



with the fermionic action 



S[c,c]=Sb[c,c]+ / cItH[c,c}, 
Jo 



(4) 



which contains the fermionic Berry term^S 

S B [C,C}=^2 d,TC a {YT)d T C a {YT) (5) 

r« JO 

that will play an important role in our considerations. 



III. SPIN-CHARGE ROTATING REFERENCE 
FRAME 



Consequently, the effective action reads: 

V(rr) V 2 (vt) 



r/3 

S[n,V,g,c,c] = V / dr 
„ Jo 



U 



U 



+ iV{rT)n(vT) + 2g(rr)f2(rr) • S(rr)] 

+ S B [c,c]+ / dTH t [c,c\. (10) 
Jo 

Since, U is the largest energy in the problem, the simple 
Hartree-Fock theory won't work. To proceed, one has 
to isolate strongly fluctuating modes generated by the 
Hubbard term according to the charge U(l) and spin 
SU(2) symmetries. 



The spin-rotational symmetry present in the Hubbard 
Hamiltonian is instrumental for obtaining proper low en- 
ergetic properties. Therefore, it is crucial to construct 
a theoretical formulation that naturally preserves this 
symmetry. In particular, one should consider the spin- 
quantization axis to be a priori arbitrary and integrate 
over all possible directions in the partition function. It 
can be achieved when the density-density product in 
Eq.(JTJ is written, following Ref. [2l], in a spin-rotational 
invariant way: 



r v 



n 2 (r-r) - [fi(rr) ■ S(rr)]' 



(6) 



where S a (rr) = §Xw c l( rr )< 5 'aa' c "'( rr ) denotes the 
vector spin operator (a = x, y, z) with a a being the Pauli 
matrices. The unit vector 

r2(rr) = [sin^(rr) cos <p(tt), sini?(rr) sin y(rr), 

cost?(rr)] (7) 

written in terms of polar angles labels varying in space- 
time spin quantization axis. The explicit spin-rotation 
invariance comes from the angular integration over f2(r-r) 
at each site and time: 



Z = [T>Sl] / [VcDc] e~ s[a '' 



where \DQ] 



sini?(rr fc )di9(rr t )dy(rr fc ) 

angular integration measure. 



(8) 



is the spin- 



A. Hubbard-Stratonovich decoupling 

The spin and charge density terms of the Hamiltonian 
in Eq.((S]) are of the fourth order in fermionic operators, 
so they must be decoupled using Hubbard-Stratonovich 
(HS) formula^ with the auxiliary fields q(tt) and iV(rr) 
respectively. The partition function can be written in the 

f orI ™23 



Z = / [Vft] / [VVVq] / [VcDc] e 



-S[n,V,g,c,c] 



(9) 



B. U(l) charge frame 

Now, we switch from the particle-number representa- 
tion to the conjugate phase representation of the elec- 
tronic degrees of freedom that is governed by the compact 
U(l) group. To this end the second quantized Hamilto- 
nian of the model is translated to the phase representa- 
tion with the help of the topologically constrained path 
integral formalism.— As a result the electrons emerge as 
composite particles consisting of of spin-carrying neutral 
fermions and topological charged bosons in a form of a 
flux tubes with the quantum phase variable dual to the 
local electron density. To this end we write the fluctuat- 
ing "imaginary chemical potential" iV(rr) as a sum of a 
static Vo(r) and periodic function 



V(rr) = Vb(r) + V(vt), 
where, using Fourier series 



oo 

V(rr) = -Y.lVirwJe^ + c.c] 



(11) 



(12) 



with oj n = 2irn/j3 (n = 0, ±1, ±2) being the (Bose) Mat- 
subara frequencies. Now, we introduce the U(l) phase 
field 0(rr) via the Faraday-type relation^ 



d^vr) = e _ i<KrT) i d e#(rr) 

dr i dr 



V(vt). (13) 



c a (rr) 




' z(rr) 




' /«(rr) " 


_ c a {rr) _ 




z(rr) 




fa(rr) 



Furthermore, by performing the local gauge transforma- 
tion to the new fermionic variables / Q (rr): 



(14) 



where the unimodular parameter |z(rr)| 2 = 1 satis- 
fies z(tt) — e i ^ rr \ we remove the imaginary term 
i Jq drV ^(rr)n(rr) for all the Fourier modes of the V(tt) 
field, except for the zero frequency. We point out here 
that a similar phase representation was developed in the 
context of Coulomb blockade in mesoscopic systems^. 
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C. SU(2) spin frame 



D. Solutions for Vo(r) and g(rr) 



Subsequent SU(2) transformation from / Q (rr) to 
h a {Yr) operators, 



R(rr) 



[/ T (rr), / ; (rr)] 

R(rr) 

with the constraint 

|C T (rr)| 2 + |Ci(rr)| 2 = l 



/i T (rr) 
h(rr) 

CtO t ) -Ci( rT 

Ci(rr) C T (rr) 



(15) 



(16) 



takes away the rotational dependence on f2(rr) in the 
spin sector. This is done by means of the Hopf map 27 



R,(rT)a z TL\rT) = a ■ fi(rr) 



(17) 



that is based on the enlargement from two-sphere S2 to 
the three-sphere S3 ~ SU(2). The unimodular constraint 
in Eq. p^j) can be resolved by using the parametrization 



Ci T (rr) = e -iM^)+x(rr)] CQS 



Ci(i 



e iW{rr)- X (vr)] ^ 



2 

tf(rr) 



(18) 



with the Eulcr angular variables (p(rr), ^(rr) and x(rr), 
respectively. Here, the extra variable x(rr) represents 
the U(l) gauge freedom of the theory as a consequence 
of S*2 — > 5*3 mapping. One can summarize Eqs p4| and 
(|15[) by the single joint gauge transformation exhibiting 
electron operator factorization 



c a (rr) = 2jW Qa '(rr)/i Q /(rr), 

a' 

W(rr) = z(rr)R(rr), 



(19) 



where U (rr) is a U(2) matrix which rotates the charge 
and spin-quantization axis at site r and time r. This 
reflects the composite nature of the interacting electron 
formed from bosonic spin and charge degrees of freedom 
given by R aa '(rr) and z(rr), respectively as well as re- 
maining fermionic part h a (rr) . Accordingly, the integra- 
tion measure over the group manifold becomes 

[Vct>vn]= V TT / # (r) fdn (r) 

{m( r)} r ^ J 
Cl(r/3)=n„ 0(r/3)=0 o (r)+27rm(r) 

x J Dn(rr) J £>0(rr),(2O) 

O(r0)=f2 0(rO)=0o(r) 

where J dfl ■ ■ ■ — Jo sm Odd J Q 27r dip . . . and 

[PJl(rr)] = Ylk d£l{vTk)- Here, m £ Z labels equivalence 
classes of homotopically connected paths 2 ^ for the U(l) 
group. 



Once can anticipate that spatial and temporal fluctu- 
ations of the fields Vo( r ) and £>( rT ) w iU be energetically 
penalized, since they are linked to the high energy scale 
set by U and decouple from the angular and phase vari- 
ables. Therefore, in order to make further progress we 
next subject the corresponding functionals to a saddle 
point analysis. The expectation value of the static (zero 
frequency) part of the fluctuating potential Vo(r) calcu- 
lated by the saddle point method to give 



IfX, 



(21) 



where /2 is the chemical potential with a Hartree shift 
originating from the saddle-point value of the static 
variable Vq(t) with = + n/i| and = 

(h a (rT)h a (rT)) . Similarly in the magnetic sector we have 



p(rr) = (-l) r A c 



(22) 



where A c = U(S z (yt) sets the magnitude for the Mott- 
charge gap. The solution delineated in Eq. (j2"2"]) corre- 
spond to the saddle point of the antifferomagnetic type 
(with staggering A c ) . Note that the notion "antiffero- 
magnetic" here does not mean an actual long-range or- 
dering - for this the angular spin-quantization variables 
f2(rr) have to be ordered as well. The mean-field param- 
eter A c has to be determined by the stationary points of 
the action, e.g., by the mean-field equations that will be 
derived more explicitly later on for our special purpose. 



IV. EFFECTIVE ACTIONS 



A. Total fermionic phase-angular action 



In the new variables the action in Eq. ([10)) assumes the 
form 

S [fl, (j>, g, h, h] = Sb [h, h] + / dTH.n,4>[p,h,h] 

Jo 



/ drA(rr) • S h (rr), 
- Jo 



where 



Sh(rr) = ^^2h a (rT)a ai h 7 (rT). 
Furthermore, 

5 o[0] = E [ dr 



(j) 2 {YT) 1 2jJL 



u 



i U 



(23) 



(24) 



(25) 



stands for the kinetic and Berry term of the U(l) phase 
field in the charge sector. The SU(2) gauge transforma- 
tion in Eq. (|15p and the fermionic Berry term generate 



5 



SU(2) potentials given by 

Rt(rr)3 r R(rr) = Rt^iL+tfil 



4 IK 



= <x ■ A(rr) 



(26) 



where 
A x {vt) 
A y (rr) 
A z {tt) 



-#(rr) sin x( rr ) — _ ( p( vr ) s i n 0{tt) cos x( rr ) 

1 . ^ 

— #(rr) cos x( rr ) + 2 ^( rr ) SU1 ^( rr ) SU1 x( rr ) 
-95(rr)cos7?(rr) + ~x(rr). 



(27) 



are the SU(2) gauge potentials. The fcrmionic sector, in 
turn, is governed by the effective Hamiltonian 



"Hq. 



(28) 



where 

Ap) 



= E(- 1 ) I ' A '=^T(rT)/iT(rr)-^(rT)/i i (rr)] 



H. 



(*) 



"* E E [^VMrV^/^rr^rV) 

(r.r') 07 

/2^/t tt (rr)/i Q (rr). (29) 



The result of the gauge transformations is that we have 
managed to cast the strongly correlated problem into 
a system of mutually non-interacting pseudo fermions, 
submerged in the bath of strongly fluctuating U(l) and 
SU(2) fields whose dynamics is governed by the en- 
ergy scale set by the Coulomb interaction U coupled to 
fermions via hopping term and with the Zeeman-type 
contribution with the massive field g(rr). 



B. Charge (phase) action 

In systems with Coulomb interactions, the phase vari- 
able dual to the charge is an important collective field. 
We start with a partition function for charge sector 



[V4>] e 



(30) 



where the charge action requires tracing over fermionic 
and angular SU(2) variables 



SI 



In / [VhVhVn]e 



-S[(p,<p,tf,h,h] 



(31) 



To proceed, it is convenient to replace the phase degrees 
of freedom by the complex field 



z (rr) 
z (rr) 



(32) 



S[c,c] 

interacting 
fermions 



z 



spin + charge rotation 

m 

S ip, 0, h, h] 

U(2) gauge fields 
+ fermions 

(5 [z,z,<;,t,h,h]) 



r 



Tr, 



L 



Tr 



SU(2) spin 
(CP') 

( S [C><1) 




S [h, h] 

rotated 
fermions 



FIG. 1: Scheme for the derivation of effective actions. The 
initial purely fermionic action of the Hubbard model is trans- 
formed into effective actions in spin, charge and fermionic 
sectors by performing the trace over selected set of variables. 



which satisfies the periodic boundary condition z (r/3) = 
z(r0). This can be done by implemented the Fadeev- 
Popov method with the Dirac delta functional resolution 
of the unityi2£ 



[VzDz] 5 

n4 z(rr) ~ 



^£i*mi 2 - 

i 

>^ (rr) l 6 h (rr) 



-i<^(rr) 



,(33) 



where we take z (rr) as continuous variable but con- 
strained (on the average) to have the unimodular value. 
We can solve the constraint by introducing the Lagrange 
multiplier A which adds the quadratic terms (in the z (rr) 
fields) to the effective action. The partition function can 
be written in form 



Z = 



/>+ioo 


\V A J 


/ —ioo 


2m _ 



(34) 



where the free energy per site T = — In Z / j3N is given 
by: 



T 



-A, 



1 

W(3 



In / \VzVz 



-S e tf\z,z] 



S c g[z,z] = E / dTdT [^*firr>8 (t - t')] 

l rr' \ 



(rr') ' 

-7 (rr, r'r')] z(tt)z(t' t'), 



where 



7 (rr, r'r') = (exp {— % [(f) (rr) — <f) (r'r')]}) 



(35) 



(36) 
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is the two-point phase correlator associated with the or- 
der parameter field, where (...) is the averaging with 
respect to the action: 



M4>} = Y, f dT 

.. Jo 



7 — 0(rr) 



[7 



i ?7 



(37) 



The action with the topological contribution, after 
Fourier transform, we write as 

5 c ft[z, z] = i ^ I (kw n ) z (kcj„) , (38) 



k.n 



where 



G zOk M=^+7 1 (<*> n ) 



(39) 



is the inverse of the propagator, while the phase correla- 
tor after Fourier transform, can be written as as: 



1 4 + °° 



:( m+ if) 2 



l-4( 



'"T[j [/ 



where 



Z = exp 



m— — oo 



(40) 



(41) 



is the partition function for the set of non-interacting 
quantum rotors. Note that the presence of the integer 
winding numbers in Eqs (|40[) and (|4ip renders the phase 
propagator periodic in the reduced chemical potential 
fi/U. The unimodular condition of the U(l) phase vari- 
ables translates into the equation 



J_v I 

Nfi f- A z0 + I' 1 K) 



(42) 



which fixes the Lagrange multiplier A 



zO- 



C. Fermionic action 

Now we turn to the effective action of pseudo-fermions 
by tracing out the gauge degrees of freedom. To this end 
we write the partition function as 



[VhVh] e 



-S[h.h] 



(43) 



The kinetic part is calculated in the cumulant expansion. 
In the first order of the expansion: 



/__« _ Jo 



dr (tg + fj,6 rT i ) 



(rr').a ' 



[/i a (rr)/i a (r'r) + ft. 



(45) 



The hopping t is renormalized by a Gutzwiller-type 
parameter— g = g c g s , where 

g a = ([Rt(rr)R(r'r)] u )={[Rt(rr)R(r'r)] 22 ) 

g c = (z(rr)z(rV)) (46) 

being a multiply of the renormalization parameters in 
charge and spin sectors. The parameters g s and g c have 
to be calculated self-consistently, they contribute to the 
band narrowing and eventually for the band collapse for 
strong correlations. Another contribution to the hopping 
amplitude comes from the the second order cumulant ex- 
pansion which generates a term of the form: 

S t (2) [ft,ft] = f dr V ^{ttt't^Ivtv't) (47) 

<rr') 

with the bond operators 

jr( rTr ' T ) = hi^hjr'T) + h l (rT)h l (r'T) ^ ^ 

The fourth-order fermionic operator terms in the action 
can be decoupled using HS transform which introduces 
additional field v. The resulting action becomes bilinear 



<S t (2) [ft, ft] = tj 



(rr')a 



/3 



dT[ft ct (rr)ft Q (r'r) + ft.c], (49) 



with the effective hopping that is proportional to the an- 
tiferromagnetic exchange constant 



t - — 



(50) 



This dispersive low-energy band, for which the band 
width is set by the exchange interaction is a clear sig- 
nature of the coupling of the quasiparticles to antiferro- 
magnetic correlations. The value of the v field can be 
fixed self-consistently using saddle-point method to give 



v = y^(ft a (rr)fea(rV)). 



(51) 



where 



S[h,h] = -In / [T>(j)T>n}e 



3 -S[v,0,i9,/i,?i] 



S B [h,h] + S? , [h, ft]+S t (2) [ft, ft] 

,(P) 



drH^[ft,ft] 



(44) 



As we will see in the following for a certain range of model 
parameters the value of v may vanish leading to the band 
collapse and the insulating state. Finally, one can write 
the resulting fermionic action in a compact Nambu form: 

S[h, h] = -^J2A h (kw n ) G^ok K) A,, (kw n ) , (52) 



7 



where the vectors (kr) are defined by 
Ah (kr) = 



hf (kw„) 
/i x (kw„) 

/if (k - 7TW n ) 

hi (k - 7ru; n ) 



(53) 



w hk 





A c 













-Ac 


A c 













-A c 








(54) 



and the inverse propagator reads 



G hOk ( w ") - 



with uj h]i = iu n - fi± 2x£k, where \ = *j/ 2 - 



D. Spin-angular action 

Since we are interested in the magnetic properties of 
the system a natural step is to obtain the effective ac- 
tion that involves the spin-directional degrees of freedom 
fi, which important fluctuations correspond to rotations. 
This can be done by integrating out fermions: 



Z 



where 



S[n] =-ln / \D(j>VhVh]e 



-S[ft] 



-S[<p,<t>,ti,h,h] 



(55) 



(56) 



generates the cumulant expansions for the low energy 
action in the form 



s[n] = s [n] + s B [n] + Sj[n]. 



1. AF exchange term 



(57) 



2. Berry term 

In general, in addition to the usual exchange term, 
the action describing antiferromagnetic spin systems is 
expected to have a topological Berry phase term 

rP 

S B [Sl] = -2V / drA(rr) • <S h (rV)), (60) 
, Jo 



where 



(S z h (rr)) = ^ 



(61) 



In terms of angular variables, the Berry term becomes 
6 f 13 

S B [Sl] = -J2 dr[^(rr)cos??(rT) +x(rr)] . (62) 
1 r Jo 

If we work in Dirac "north pole" gauge x( rr ) = ~ ( fi( rT ) 
one recovers the familiar form 

a rP 

S B [fl] = -V/ dT<p(rr)[l - cos-&(tt)]. (63) 
1 ^ Jo 



Here, the integral of the first term in Eq. (j62|) has a 
simple geometrical interpretation as it is equal to a solid 
angle swept by a unit vector f2(#, ip) during its motion. 
The extra phase factor coming from the Berry phase, re- 
quires some little extra care, since it will induce quantum 
mechanical phase interference between configurations. In 
regard to the non-perturbative effects, we realized the 
presence of an additional parameter with the topological 
angle or so-called theta term that is related to the Mott 
gap. In the large-Z7 limit one has A c — > £7/2, so that 
— > i relevant for the half-integer spin. However, for ar- 
bitrary U the theta term will be different from that value, 
which, as we show, will be instrumental for destruction of 
the antifferomagnetic order away from the spin-localized 
U —* oo limit. 



The part of the action that involves the spin stiffnesses 
is given by 

Sj[Sl} = [ dTfl(r T ) -flir'T) (58) 

<rr') J ° 

with the AF-exchange coefficient 

r. s 4t 2 , , 2 At 2 /2A C \ 2 , , 

J(A C ) = — {n^nif^ — f-^J . (59) 



From the Eq. 



it is evident that for U 



oo one 



has J(A C ) ~ ^j- since _> \ m this limit. In general 
the AF-exchange parameter persists as long as the charge 
gap A c exists. However, J(A C ) diminishes rapidly in the 
U/t — ► weak coupling limit. 



3. Kinetic energy term for spin 

In analogy to the charge U(l) field the SU(2) spin sys- 
tem exhibit emergent dynamics. Integration of fermions 
will generate the kinetic term for the SU(2) rotors 



rP 

S [n] = 2 V / dr 
- Jo 



XT £ A a {vr)A a {vr) + 



a—x,y 

+ X lA z {vt)A z {vt)], 
which can be written in a more compact form: 



(64) 



rP r 

S [Q] = 2]T / dr xt V] A (rr) • A (rr) + 

r J L 

+ (xl-Xt)A z {tt)A*{tt)], (65) 
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where: 



Consequently, the spin-angular action transforms into 



XT = (SI{tt)SI{tt)) = {SI(tt)SI{tt)), 

Xl = (rr) (rr)) . (66) 

The transverse susceptibility behaves in weak and strong 
coupling limit as follows^ 



XT 



8J 

_Li/T t>>a 



(67) 



Thus, in the large-!/ limit the spin kinetic part vanishes 
and the Hubbard model maps, as expected, onto a quan- 
tum spin- 1/2 Heisenberg model with near neighbor an- 
tiferromagnetic exchange integral J = 4t 2 /U. Here, the 
superexchange interaction J is determined by the virtual 
hopping of an electron of a given spin to an adjacent site 
containing an electron with an opposite spin. Thus the 
dynamics of J involves virtual excitations above the Mott 
gap which is set by U, and the effective interaction is es- 
sentially instantaneous. However, the retarded contribu- 
tion occurs on an energy scale which is small compared 
to the bare bandwidth and the onsite Coulomb interac- 
tion. For the correlated systems, the relative weight of 
the retarded and nonretarded contributions to the effec- 
tive interaction remains an open question. 



4- CP representation 

Now, we use a compact matrix notation adapted to the 
SU(2)-invariant character of the Hamiltonian and the ef- 
fective action including a consistent scheme of coherent 
states within a functional-integral formulation, ft the 
CP 1 representation the spin-quantization axis, can be 
conveniently written as 



n (rr) = ]T ( a (rr) <r aa , ( a , (rr) . (68) 



Therefore, all the terms in the spin action can be ex- 
pressed as functions of £ Q (rr), ( a (rr) variables instead 
of angular variables, which are more complicated to be 
handled. The spin-kinetic and Berry phase term now 
assume the simpler form 

A (it) • A(rr) = — - ?9 2 (rr) + ip 2 (rr) 

+ x 2 ( rr ) + 2</j(rr)x(rr) cos i?(rr)l 
= -C (rr) • C (rr) , 

A z (rr) = -<p{tt) costf(rr) + -x(rr) 



l r 
2 



S [C, C] = V [ dr \2 X tC (rr) • £ (rr) 
[C(rr)-C(rr)-C(rr)-C(rr) 



XL - XT 



- 6 (-l) r [C (rr) • C (rr) - C (rr) • C (rr)] } 
-J^2 / dTA(Trr f T)A (rrr'r) (70) 

with the bond operators: 

A (rrr'r) A (rrr'r) = ~Sl (rr) • ft (r'r) + i 
^ (rrrV)= CTMa(rV ) ^q(rr ) C T (r'r ) (n) 

relevant for the bosonic representation of an 
antiferromagnet ^ 



5. Canonical transformation of CP 1 variables 

In order to achieve a consistent representation of the 
underlying antiferromagnetic structure, it is unavoid- 
able to explicitly split the degrees of freedom accord- 
ing to their location on sublattice A or B. Since the 
lattice is bipartite allowing one to make the unitary 
transformation^ 



C T (rr) - -Cx(rr) 
Cx(rr) - C T (rr) 



(72) 



for sites on one sublattice, so that the antiferromagnetic 
bond operator becomes 



^4(rrr'r) -> A' (rrr'r) = 



Ca(rT)Ca(r'r) 



V2 



(73) 



This canonical transformation preserves the constraint 
in Eq. (|16[) . Biquadratic (four- variable) terms in the 
Lagrangian cannot be readily integrated in the path in- 
tegral. Introducing a complex variable for each bond 
that depends on "imaginary time" Q(rrr'r) we decouple 
the four-variable terms ^4'(rrr'r)^l'(rrr'r). In a simi- 
lar manner by introducing a local real field a (rr) , we 
can decouple second term in the r.h.s. in the Eq. ([70)) . 
To handle the unimodularity condition one introduces a 
Lagrange multiplier Af (r). to treat the variables ( a (rr), 
C a (rr) as unconstrained bosonic fields. Consequently the 
effective Hamiltonian becomes 



8 rT > 



C(rr)-C(rr)-C(rr)-C(rr) 



(69) 



(rr') 

+ QC-C + QC-C + kC-CS rr ,] (74) 
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where 



A. Charge sector 



XT - XL 



(75) 



The saddle-point values of the Q, a and fields are given 
by 



a sp (rr) = X <C • C - C • C> =0 

? sp (rrr'r) = ~{C(rr) ■ C(rV)) 
1 = (C(rr) • C(rr)) 



(76) 



and by assuming the uniform solutions Q sp (rrr'T) = Q, 
a sp (ttt't) = a and A^ sp (r) = A^ we obtain for the 
Hamiltonian in the spin-bosonic sector 



H [C. C] = ^ E A Cff (kw n ) G^ k K) A fCT (kw n ) 



knu 



with 



(77) 



With the charge sector Green's function from Eq. ([39| 
one can write the constraint for Lagrange multiplier A z : 



cothf^-cothf^ 

P „ A.I f (M\ 2 i 25A 



where 



(80) 



2<5A 2 



(81) 



with f{2(j,/U) = frac(2/i/£7- 1/2) - 1/2, where 
frac (x) = x — [x] is a fractional part of x, since for tem- 
peratures much lower than U (/? -C f ) the summation 
over winding numbers leads to periodic dependence of 
the model on the chemical potential. 



B. Fermionic sector 



and 



A C<T (kw n ) 



_Ct (k,w n ) 
Co- (-k,-w n ) 
_C CT (k-7r,w„) 

C CT (~k + 7T, -u„) 



(78) 



^ + A C 2Q£ k -2i0w„ 

-2i9uj n 

2i0u;„ -2Q£ k ^ + A c 



+ X c 2i6»w„ 
+ \ — 2Q£k 



(79) 

and £ s = 1/ (2^r) sets the kinetic energy scale for the 
SU(2) rotors. 



The fermionic sector introduces the Mott-charge gap 
A c and the v field, which renormalizes hopping in the 
second order of cumulant expansion: 

k 

k 

where njr (£") is the Fermi distribution, 



V. SELF-CONSISTENT EQUATIONS 

Procedure of decoupling of the Hubbard Hamiltonian 
in Eq. (fTJ) introduces numerous decoupling fields, which 
values are fixed within saddle-point approximation. Also, 
the bosonic degrees of freedom in charge and spin sectors 
can condense leading to superconducting and magnetic 
ordering (although in the present paper we only consider 
magnetic ordering within the spin sector). Since, the 
ordering is described within quantum rotor model, it in- 
troduces additional constraints in each sectors for order 
parameters or Lagrange multipliers. Together, it creates 
a set of non-linear self-consistent equations, which allow 
for calculation of the effective variables of the present 
theory. 



E± = -/l ± E k , with £ k = J A* + ( ) ! ^ ) 



and the lattice structure factor 



Ck = cos k x + cos k y . 



(84) 



The presence of A c gives rise to a Fermi surface insta- 
bility as first suggested by Slater^ The basic principle 
behind an antiferromagnetic Slater insulator can be ex- 
plained most easily by electrons living on a bipartite lat- 
tice, that could be separated into two inter-penetrating 
sublattices (let say, A and B) such that the nearest neigh- 
bor of any site are members of the opposite sublattice. 
In the corresponding band structure picture the lattice 
unit cell is doubled and the first Brillouin zone is cut in 
half. 
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C. Spin sector 

In the Hubbard model the quantum-mechanical ob- 
jects are not local spins but mobile electrons such that 
we have to expect that the analysis of the ground-state 
phase diagram as a function of the interaction strength is 
even more difficult than for the Heisenberg model. So far 
in the above discussion we did not consider the possibility 
of an ordering of magnetic moments. Thus, the concept 
of the Slater insulator has to be supplemented by that of 
the Mott-Heisenberg insulator displays long-range order. 
Since correlations are absent in a Hartree-Fock descrip- 
tion the moments order at the very same temperature 
they are formed. In contrast, the moments are already 
present in the Mott-Heisenberg insulating state and re- 
main in the paramagnetic phase. Therefore the order- 
ing of the pre- formed moments, the important signature 
of electron correlations, provide a clear distinction be- 
tween the ideas of Slater (self-consistent single-electron 
theory) and Mott (many-electron correlations). There- 
fore, a nonzero value of A c does not imply the existence 
of AF long-range order. For this the angular degrees 
of freedom ri(rr) have also to be ordered, whose low- 
lying excitations are in the form of spin waves. In the 
CP 1 representation (where the Neel field is represented 
by two Schwinger bosons) Bose-Einstein condensation of 
the Schwinger bosons at zero temperature signals the ap- 
pearance of AF long-range order. The AF order param- 
eter in terms of the original fermion operators is defined 
as 

mAF = 5>l) r (^(rr)> = 

r 

= ^(-l) r (n(rr)).(S,(rr)). (85) 

r 

Owing the fact that (S%(tt)) = (— l) r A c 5 Qj2 we obtain 
m AF = A c ^(fF(rr)) 

r 

= A c £ [(C T (rr)C T (rr)) - (Ci(rr)Ci ("")>] (-86) 

r 

Furthermore, the order parameter for the CP 1 "boson 
condensate" 

(Ca(ko; n )) = (Ca(kw„)) 



where the equation fixing the order parameter TOq reads 



j3N 

—^-moSo,u n Si, a (<Sk,o + £k,7r) ■ (87) 

This yields a macroscopic contribution (i.e., order one) 
to the staggered magnetization and represents a macro- 
scopic contribution to the CP 1 bosons density too, of the 
a =| bosons at the mode with k = 0, uo n = thus giving 

m AF = J^Yl [(Ct( kw ™)CT( ka, n)) 

-(Ci(ko; n )Cx(ka; n ))] = ^m 2 , (88) 



1 — TOn 



J2 G ck K 



(3N 



1 \ - f coth (lj- )] + coth [(3E\ 



k V^k. 



coth [0E- k (w£)] + coth [(3E\ 



sk V^k , 



(89) 



Also, decoupling of the bond operators in the kinetic term 
of the spin action in Eq. ([70)) leads to additional field Q, 
which value is determined from the equation: 



2J 

z(3N 



7n 



zN 



+ 



coth [pE- k (cj k )] + coth [0E+ K)] 



(90) 



coth ( w +)] + coth [0E+ ( W +)] 



where 



and 



cj± = a c ± 2Q&. 



(91) 



(92) 



From the Eq. (|89|) it follows that the magnetic order at fi- 
nite temperatures is excluded in two dimensions in agree- 
ment with Mermin- Wagner theorem. 33 



VI. SINGLE-PARTICLE SPECTRAL 
FUNCTIONS 

We turn our attention now to the spectral function de- 
fined in terms of the electron Green's function. Within 
our construction, it is possible to write the electron 
Green's function as a product of U(l) phase, SU(2) 
spin (in CP 1 representation) and pseudo- fermion Green's 
functions: 

G z (rr,rV) = - (z (r,r) z (r',r')> , 

Gf(vT,v'r') = -<C«(r,T)Ca'(r , y)>, 

GrV,rV) = -(h a (r : r)h al (r\r')), (93) 

while the full Green's function of the system is the prod- 
uct: 



G aa (rr'r) = - £ (z r z T ,) {R aP (rr) R\ a (r'r)) 

x {hp (rr) h-y (r'r)) . (94) 
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FIG. 2: (Color online) Evolution of the density of states of 
the model for various interactions U/t from strong to weak- 
coupling limit. 





U/t=4 








U/t=1.49 



3 

>■ — ' 



U/t=1.22 




FIG. 3: (Color online) Evolution of the pseudo-fermionic den- 
sity of states [see, Eq. (|B9[) ] of the model for various interac- 
tions U/t from strong to weak-coupling limit. 



The problem of calculating the spectral line shapes now 
becomes one of calculating the convolution of the Green's 
functions in Eq. (|93p . Since, in the antiferromagnetic 
phase we allow for ordering in the spin sector, the aver- 
ages over spin variables can be non-zero since CP 1 Bose 
condensation signals the appearance of AF long-range 
order 

(Ca (k,w„)) = (Cq (k,w n )) 

f3N 

-^-rn S 0}UJn S^ a (<5 k ,o + ^k,7r) ,(95) 

where the order parameter mo measures the fraction 
of the condensed CP 1 bosons. Consequently, the spin 
Green's function can be split into two contributions ac- 
cording to 

Gfk K) = - (Ca (kWn) Ca (kw„)) 

= - [(Ca (00) Ca (00)) + (Ca (ttO) Ca (w0))] 

} (C Q (kw„) C« (kw„)) . (96) 



Because, in this case the average is periodic with respect 
to k: 

(Ca (kw„) C~ a (k£J„)) = (Ca (k + 7T£J„) Ca (k + 7TtJ„)) , 

(97) 

the spin sector Green's function reads: 

G"k k) = -/3Wm 2 <Mko<L„o 

+ (l-2tf QT S ko <L„o)G<k (98) 

Substituting this result into Eq. (fM| one can calculate 
spectral density of the system (for details, see Appendix 
0: 

( W ) = mgASa, (w) 
x [ra B (-W + w') + n F (a/)] 

(J 

and the density of states: 

PT (w) = mlp«% (to) + 2p?% (w) . 



(99) 



(100) 
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w/t 

FIG. 4: (Color online) Spectral weight relocation induced by 
change of U /t interaction at zero temperature. The weight of 
coherence peak (equal to mj) present for U jt — 1.2 spectrum 
is relocated to the higher frequencies for U/t = 1.178, so the 
norm (spectrum integrand over frequency) is preserved and 
equal to 1. 

The single electron density of states contains then two 
terms. The first generates the the coherence peak asso- 
ciated with the log-range AF order, which is a product 
of the condensate density and the pseudo-fermion spec- 
tral function. To the extent that the peak and the back- 
ground are distinguishable objects, see Fig. 2, the weight 
under this quasi-particle peak should be proportional to 
the condensate density of the CP 1 bosons. Additionally, 
we have calculated the density of states corresponding 
only to fermionic propagator G^ Q (ki/ n ) [see, Eq. (|B9|) ] to 
illustrate the correspondence between Hartree-Fock ap- 
proach and our method. The outcome is presented in 
Fig. 3: the charge gap is a monotonic function of in- 
teraction strength U, which is the feature of Hartree- 
Fock approaches. 34 Upon crossing the antiferromagnetic 
phase boundary, one observes a remarkable feature in 
the spectral density plot. Figure 4 shows how the spec- 
trum evolves upon entering the ordered magnetic state 
as a function of the Coulomb interaction. In a conven- 
tional Fermi liquid system, the spectrum would contain 
a narrow peak, which is the signature of a well-defined 
quasiparticle. In the AF state exactly such a peak is 
seen. As Fig. 4 shows, the quasiparticle peak disap- 
pears upon leaving of the ordered state. At the same 
time the gap is preserved, yet the narrow peak is gone. 
As the peak disappears the spectral weight is transferred 
to the incoherent excitation background. Evidently, the 
quasiparticle owes its existence to the of the condensa- 
tion of the CP 1 bosons in the antiferromagnetic state, 
and not the energy gap. This exhibits a great similar- 
ity to that, which is seen in the ARPES spectra in the 
underdoped samples, however in this case the conden- 




-3 -1.5 1.5 3 

w/t 



FIG. 5: (Color online) Evolution of density of states as a 
function of the temperature (as indicated on the plot) for of 
U/t = 1.25. 

sate density refers to the superconductor. The electron 
spectral function is very broad above the superconduct- 
ing transition, but a sharp quasiparticle peak develops 
at the lowest binding energies, followed by a dip and a 
broader hump, giving rise to the so-called peak-dip-hump 
structure^ 5 , which is very similar to our results depicted 
in Fig. 4, where the spectral density and density of states 
are sums of a coherent part consisting of convolved func- 
tions, proportional to the order parameter and incoher- 
ent - being a convolution of charge, pseudo-fermion and 
spin functions. It is interesting to note that the evolution 
of a Mott-Hubbard insulator into a correlated metal has 
been examined in the two-dimensional Hubbard model by 
using the cellular dynamical mean-field theory, 36 which 
incorporates short-range spatial correlations. At half fill- 
ing these correlations create additional bands due to the 
ordered antiferromagnetic states that bear similarity to 
our results. As far as the comparison with earlier works 
based on perturbation theory 4 ' 5 is concerned, the main 
effect of correlations in the form of a transfer of the spec- 
tral weight to the high energies is reproduced. The evolu- 
tion of the spectral density as a function of temperature 
is depicted in Fig. 5. At finite temperature there is no 
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AF ordering according to the Mermin- Wagner theorem 
and consequently no coherence peak, but one observes 
the gap filling as the temperature increases. 



VII. CONCLUSIONS 

In this paper we have presented a method of calcu- 
lation of spectral densities for strongly correlated sys- 
tems in terms of a collective phase variable, the rotating 
quantization axis and the fermionic degrees of freedom. 
In systems with strong Coulomb interactions, the phase 
variable dual to the local charge is an important collec- 
tive field. A theory of the Hubbard model involving free 
fermionic degrees of freedom self consistently coupled to 
a quantum U(l) and SU(2) rotor model has been devel- 
oped. The most interesting aspect of our approach lies 
however in the possibility of going beyond purely local 
mean-field description by incorporating the effect of spa- 
tial correlations and in particular the influence of the 
ordered states on the spectral properties of the system. 
This is very similar to the method implemented in Refs. 
Il7lll8l . where the single-particle properties are obtained 
by writing the fermion field in terms of a Schwinger bo- 
son and a pseudofermion whose spin is quantized along 
the fluctuating Neel field. However, in the above works 
the charge sector was treated on the mean-field level only 
and the lattice was approximated by the continuous long- 
wave limit. In our approach, the inclusion of the antifer- 
romagnetically ordered phase was done by resorting to 
the saddle-point analysis of the bosonic and fermionic ef- 
fective actions, however the general architecture of the 
method is not resting on this assumption. The method 
is suitable for general value of the correlation energy U, 
which is in contrast to the TPSC metho d 14 ! 15 valid in a 
weak-coupling limit and slave boson (fermion) approach, 
where the normalization of the spectral function is also 
violated. 37 In the large-?/ limit, the theory is controlled 
by the parameter J — At 2 /[/, in agreement with previous 
approaches. Regarding the critical behavior, in the vicin- 
ity of the phase transition, our model based on the spher- 
ical approach will be equivalent to that of n- vector model 
in the limit n — > oo, which is the same as in the TPSC 
approach. We investigated the transformation properties 
of Hubbard Hamiltonian leading to a symmetry adapted 
formulation, which explicitly exhibits the SU(2) and U(l) 
invariance included in the formalism. In this picture, the 
collective bosonic modes that represent charge and spin 
play an important role since the are related to the under- 
lying symmetries of the system and the ordered states. 
The single-particle properties are obtained by writing the 
original fermion field in terms of a U(l) phase field related 
to the charge, CP 1 bosons that parametrize the variable 
quantization axis related to the rotational symmetry and 
a pseudo-fcrmion. This decomposition allows us to write 
down the fermion Green's function by the product in real 
space of the phase, Schwinger boson propagator and the 
remaining fermionic propagator. Because spatial correla- 



tions are now included, we find important modifications 
of the electronic picture due to the formation of the or- 
dered magnetic states. We have shown that the single 
electron density of states, consists of two pieces. The first 
generates the the peak which is a product of the conden- 
sate density of CP 1 bosons that represents the antiferro- 
magnetic order and the and the pseudo-fermion spectral 
function. We found that this feature is an analog to the 
situation present in the normal state ARPES spectra in 
the underdoped samples, where the the superconducting 
condensate produces similar behavior. 

Finally, it would be interesting to extend our calcu- 
lation to the more interesting doped system, which will 
require additional numerical effort. 
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APPENDIX A: CORRELATION FUNCTIONS 

1. Charge sector 

The Green's function in the charge sector reads: 

1 



G z (u„) 



A* + f + £k-^/(^)]' 



(Al) 



Here, function / (2///C7) = frac(2/i/C7) — 1, where 
frac (x) = x — [x] replaces summation over winding num- 
ber in Eq. pTJ|) , which is valid for temperatures (3 < [/ '. 



2. Fermionic sector 

In the fermionic sector, the Nambu notation of the 
fermionic action in Eq. (|52p with vectors: 



A(kw„) = [/ltk,ft-lk,/lTk-7r,ft-lk-7r] k) 

leads to the Green's function matrix: 



(A2) 



u hk 





r hk 








/-til 

^hk 





r hk 


r hk 





r ll 

^hk-TT 








ipll 

r hk 








K) (A3) 



with normal and anomalous Green's functions 

A* + 2 x£k - iv n 



Fr (kl/n) 



A? + «2_ (Ai _ Wn ) 2 ' 



(-lfA c 



A2 + 4 X 2 ^-0u-^ 



(A4) 
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3. Spin sector 

In the spin sector, the spin action in Eq. (|77|) with 
vector 

A(k,CJ n ) = [C CT k,Ca-k,Cak-7r,Ccr-k+7r] k) (A5) 

leads to spin Green's function matrix: 



which elements read 
G^k (w n ) 

-F?k K) 
-F?k (w„) 



G^k, ^Ck, 


_G^k, 




^Ckj _^Ck) 


-F^Ck— n-, 


G(Jk— 7T 


G^k, ^Ck-7r, 


G(Jk— 7T, 


-FCk— 7T 


. -Fck; G^k-7r 


-^Ck— !r 


G|Jk— 77 . 



(A6) 



- [7_ (k,w„) +7+ (k,w„)] 

- [7_ (k,w n ) - 7+ (k,w„)] 



7 g (k,w n )+7 g (k,w„) 
7 Q (k,w n )- 7 ?(k,w n ) 



(A7) 



where 



7± (k, Ul r , 



^ ± 2i6u n + A 



y 2 - ± 2i6uj n + A^ 
2Q& 



^ ± 2i(9w„ + A 



(A8) 



APPENDIX B: SPECTRAL DENSITIES 

The Green's function of the system is the combination 
of Green's functions of charge, spin and fermionic sectors: 

Gr'(rr,rV) = S aa ,G z (rr, rV) 

+ Gf (vr,v'r')Gf (rr,rV)],(Bl) 
where its Fourier transform: 



x G QQ (r-r',r-r'). 



(B2) 



The spectral density is the imaginary part of the single- 
particle Green's function and therefore contains full in- 
formation about the temporal and spatial evolution of a 
single electron or a single hole in the interacting many- 
electron system. The spectral density is defined for 
fcrmions as follows 



G a x l K) 



+ °° du, A$° k H 

^ 2ir iv n - lo 



(B3) 



with X = c, h for full system and fermionic sector, re- 
spectively. Similarly, for bosonic sector: 



Gjfk (w n ) 



dw A X k (u>) 
2tt iu> n — IjJ 



(B4) 



where X = z, £ for charge and spin part. The full spectral 
function of the system expressed in Fourier variables is a 
double convolution of three elementary spectral functions 
of charge, spin and fermionic sectors, which written in 
terms of the real frequencies are 

-y 



Aaa 



x [n B (-u)') + n F (oj - u)')] 

U ) = N^J ^ t^Cq^ ' .<>k- (1 

x [n B (-UJ 1 ) + n F (ui - ui')] 



(B5) 



Introducing a density of states being a local (k- 
intcgrated) spectral density defined as 



2ttN 

k 

one obtains convolution expressions 

e + OG 



(B6) 



J — oo 

x [n B (-a/) + uf (o> - w')] 

/+oo 
duj'p^ ( w ')p£( w - w ') 
-oo 

x [n_B (— to') + n F (w — uj')] ■ 
The density of states in the charge sector reads: 

Pz ( W ) = *("- Sc-) 

while in the fermionic sector 
f/ 2 |w-/2| 



(7 



(B7) 



(B8) 



-e 



U 4 v 2 (B+ - B 
x [P2d{B-)+ P2d{B+)] 



with 



B±=±- 



UyJ( u-ji) -Ag 
2i 2 v 



(B9) 



(BIO) 



In a non-dispersive case (t = or u = 0), the fermionic 
density of states 

P r H = is (u - p + A c ) + 1 -8{uj- p - A c ) . (Bll) 
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In the spin sector, 



4Q 



sgn 



P2D 



A c - 29lo 



2Q 



sgn 



P2£> — 



t- - A c + 26»w 



2Q 



where 



P2D (x) 



1 



K 



1-f 



(B12) 



(B13) 



energy w to a many-body system. Therefore, correctly 
calculated spectral function of the system should meet 
the normalization condition. We can verify this is indeed 
the case for the scheme presented in the present work. 
Since, (uj) is given by the Eq. (|99j) . its norm reads: 



dco 
2^ 



^A^{ U ) + ^A%Z Ck (uj) 



(C2) 



is the density of states for the square lattice and K(x) 
stands for the complete elliptic integral of the first kind. 38 



APPENDIX C: NORMALIZATION CONDITION 



Considering the first of two terms ^4"^ k (uj) and using 
the expression for the convolved charge-fermionic spec- 
tral density from Eq. (|B5[). one obtains: 



Among the general properties of the spectral function 
there are several sum rules. A fundamental one is 



2^ 



A% (uj) = 1. 



(CI) 



which reminds us that (uj) describes the probability 
of removing/ adding an electron with momentum k and 



Transforming the integration variable u> — > uj — uj' we 
obtain: 



r°° rf ^« a/ , 1 r f ^ a , is , f duj A« a , \ 



r 



Since, the charge spectral function is antisymmetric mation rule: 
A Z q (uj) = —A zcl (— w), its integral over frequencies van- 
ishes: 



+ OC 



2^ 



A zcl (w) = 0. 



(C5) 



On the other hand, the integrand of fermionic spectral 
density is normalized 



+ OC 



2^ 



Ag£_ q (w) = l. 



The remaining factor: 
; duJ ' 



q'#0 



2?r 



A ZCL , (uj')n B {-uj') = 1 



(C6) 



(C7) 



is simply a charge sector constraint, which can be checked 
by substituting Eq. (|B4|) to Eq. (|5D|) and using a sum- 



1 y- 



o i0~ 



(3 ^ iuj n - uj e-^ - 1 
Finally, the first term in the Eq. (|C2[) reads 

J~oo 27T 



(C4) 



n B (-uj) . (C8) 



(C9) 



Similarly, one have to treat the second term in the Eq. 
(IC21) 



x^ShU, (w - w') [n fl (-«/) + n F (uj - uj')} . (CIO) 
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Once again A"^ (u>) — —A^ (— a;), which leads to the 
vanishing of the integral 

/ + OO 1 
— ^ C qH=0. (Cll) 
-oo 27r 

The remaining part is the spin sector constraint 

^EE r d ^ { J)n B ^') = l-ml 

(C12) 

It means that the second term in the Eq. (|C2[) is equal 
to 

/ + OO J 
^ — A^ fk (a;) = l-mg. (C13) 



Substituting the results to the Eq. (|C2[) one can see that 
the normalization condition from the Eq. (|Cip is always 
fulfilled. The norm of the full spectral density directly 
depends on the constraints in the charge and spin bosonic 
sectors. Therefore, careful solution of the self-consistent 
equations in Sec. |V]is of primary importance to obtain 
physically reasonable results. 
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